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Please  read  the  following  problem  formulation  and  select  among 


the  options  listed  below. 

A  warehouse  contains  N  batteries  in  storage.  Batteries  fail 
(silently)  according  to  life  distribution  F(t)  =  1  -  exp(-Xt),  X  known. 

At  time  tj  we  select  a  random  sample  of  size  n^  batteries  and  find 
that  fj  are  failed  (and  discarded  immediately),  while  n^  -  f^  are 
functioning  and  returned  to  the  warehouse.  At  time  t^  >  tj  we  are  to 
select  a  second  random  sample  of  n^  batteries.  How  shall  we  select  our 
second  sample:  more  specifically,  how  many  as  yet  untested  batteries 
shall  we  select  and  how  many  batteries  tested  at  time  t^  shall  we  select? 
(Of  course,  the  total  number  selected  must  equal  n^.) 

1.  Select  as  many  from  among  the  untested  batteries  as  possible. 


If  the  sample  size  n2  is  bigger  than  the  number  of  such  batteries, 
choose  the  remainder  from  among  the  previously  (at  time  t^)  tested 
batteries. 


2.  Select  as  many  from  among  the  previously  tested  batteries  as 
possible.  If  n^  >  n^  -  f^,  to  complete  the  sample,  choose  the 

remainder  from  among  the  untested  batteries. 

_  _  .  -At?  .  _ -A(t0  -  t,).  If  p-  is  closer  to 

3.  Compute  pt  =  e  z  and  p2  =  e  *•  2  V  v2 

1/2  than  is  p^,  use  the  sampling  procedure  described  in  2.  If 
Pj  is  closer  to  1/2  than  is  p2<  use  the  sampling  procedure  des¬ 
cribed  in  1. 

4.  Dismiss  the  problem  formulation  as  unrealistic,  since  who  knows 

X?  Set  up  a  Bayesian  model  and  proceed  as  in  4 '  below. 

You  guessed  it  -  all  four  choices  are  correct.  That  is,  there  are 

three  reasonable  optimality  criteria  such  that  for  a  given  criterion,  the 
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optimal  procedure  is  described  in  one  of  the  first  three  statements 
above.  The  Bayesian  procedure  proposed  in  statement  4  is  appropriate 
for  X  unknown. 

Let  us  take  the  easy  cases  first.  Suppose  our  goal  is  to  maxi¬ 
mize  the  number  of  batteries  that  we  know  are  functioning  at  time  t^. 

Then  it  is  easy  to  see  that  we  should  follow  Procedure  2.  The  proba¬ 
bility  that  a  battery,  found  to  be  functioning  at  time  t^,  will  be  found 
upon  inspection  at  time  t^  to  be  functioning  is  e  **-/ e  ***  3  e  ^ 

Cn  the  other  hand,  a  previously  uninspected  battery  has  probability  e  **2 
of  being  in  the  functioning  stai;?  upon  inspection  at  time  t^.  Since 

it  follows  that  the  optimal  strategy  to  maximize 
the  number  of  batteries  we  know  to  be  functioning  at  time  t2  is  as  des¬ 
cribed  in  Procedure  2. 

Suppose,  on  the  other  hand,  our  goal  is  to  weed  out  as  many  defective 

batteries  as  possible  by  inspection.  Now,  the  optimal  plan  is  described 

in  Procedure  1.  This  follows  from  the  fact  that  the  probability  that  a 

battery  found  to  be  functioning  at  time  t^,  will  be  found  upon  inspection 

at  time  t^  to  be  in  the  failed  state  is  (e  -  e  **2)/©  ^!  *  1  -  e  *^*1  t2^. 
On  the  other  hand,  a  previously  uninspected  battery  has  probability  1  -  e  **2 

of  being  in  the  failed  state  upon  inspection  at  time  t2*  Since 

1  -  e"**^  >  1  -  e  ^*2  "  it  follows  that  the  optimal  strategy  to 

weed  out  the  maximum  number  of  defective  batteries  is  as  described  in 

Procedure  1. 

A  third  possible  goal  is  to  estimate  as  precisely  as  possible  -  precision 
is  defined  by  the  squared  error  loss  function  -  the  number  X  of  functioning 
batteries  at  time  t ^  when  X  is  known.  The  problem  is  to  select  a  sample 


of  fixed  size,  from  among  (a)  the  batteries  already  inspected  at  time 


tj  and  found  to  be  still  alive  and  (b)  the  batteries  not  yet  inspected 
at  time  tj. 

Note  that  the  (Bernoulli)  variable  used  to  indicate  the  state 
(functioning  or  not)  of  a  battery  selected  from  among  (a)  has  variance 
e"^(t2  ”  *l'(l  -  e-X(t2  "  *1^),  while  the  •‘viable  used  to  indicate  the 
state  of  a  battery  selected  from  among  (b)  has  variance  e  Xt2(l  -  e  Xt2). 
Since  we  wish  to  obtain  an  estimate  of  the  number  X  with  smallest  possible 
risk,  we  simply  follow  procedure  3  above. 

4'  Suppose  that  X  is  unknown  and  we  choose  a  distribution  to  repre¬ 
sent  our  prior  opinion  about  X.  In  particular,  suppose  we  choose  a  Gamma 
distribution  a  priori.  That  is, 


ft0!®*1  -AX 

n(X)  *  rfe)  e  (1) 

is  the  prior  density  for  X.  Having  observed  a  sample  of  size  n^  at  time 
t,  containing  fj  failed  batteries  and  n^  -  f^  functioning  batteries,  we 
compute  the  posterior  mean  of  the  function: 

Pj(l  -  p)  -  P2(l  *  P2)  *  e~Xt2(l  -  e‘Xt2)  -  e  X^2  tP(l  -  e  X^2  *1^)  * 

e-»t2  _  e-2X,.  ->'VV  ♦  (2) 

where  p^  and  p^  are  defined  in  statement  3  above.  If  the  result  obtained 
is  greater  than  zero,  we  select  one  battery  from  among  the  N  -  n^  completely 
untested  batteries  and  test  it  at  time  tjl  if  the  result  obtained  is  less 
than  or  equal  to  zero,  we  test  at  t2  a  battery  chosen  from  among  the 
n^  -  fj  batteries  that  survived  the  test  at  time  t^.  We  incorporate  this 
new  result  in  the  original  sample  and  obtain  an  adjusted  likelihood.  With 
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this  likelihood  and  the  original  prior  we  again  compute  the  posterior 
mean  of  (2)  in  order  to  decide  from  which  set  (a)  or  (b)  the  next 
battery  will  be  selected.  This  procedure  will  be  followed  sequentially 
to  complete  the  sample  of  size  n2  or  until  either  sets  (a)  or  (b)  are 
exhausted;  in  this  case  complete  the  sample  (of  size  n2)  with  items 
chosen  from  the  remainder  set  (either  (a)  or  (b)). 

To  specify  the  computations  explicitly  we  list  the  following  notation: 
nj  ■  number  of  batteries  tested  at  time  tj. 
fj  ■  number  of  batteries  (among  these  n^)  failed  at  time  t^. 
nj2  ■  number  of  batteries  tested  at  time  tj  and  at  time  t2  • 
f j2  »  number  of  batteries  failed  at  time  t2  from  among  these  n^2 . 
n2  *  number  of  batteries  tested  at  time  t2  . 

f^  ■  number  of  failures  among  the  n2  -  n^  batteries  tested  only  at  time 
t2>  (Note  that  the  f ^  batteries  are  not  part  of  the  f 2  batteries.) 

For  reasons  of  simplification  we  relax  the  above  notation  in  the 
following  formulas.  The  quantities  n2,  nj2,  f2.  and  f ^  are  viewed  below 
■;3  chan^i: o  stvR^ily  r.t  c'ch  st.-p  of  the  sequential  procedure. 

The  general  likelihood  is  given  by: 

L  •  (e‘Xtl)nl_n12‘fl  (i  .  e‘Xtl)fl(e‘Xt2)n2*f2‘f12  (1  -  e'Xt2)f2 
*  (e'Ul  -  e'Xt2)f12  (3) 

We  emphasize  the  fact  that  this  likelihood  is  computed  at  each  step  with  the 
current  values  of  n2,  n12,  f2»  and  f 12» 

From  the  likelihood  (3)  and  the  prior  (1),  we  may  compute  the  Bayes 
estimator  for  e"X^c2t2"cltP ,  as  follows: 
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Note  that  the  final  Bayes  estimator  for  e  **2  is  computed  by  taking 
c2  ■  1  and  Cj  ■  0  in  (4) . 

He  emphasize  that  the  Ba;'es  sequential  procedure  described  above  may 

•  ** 

be  only  an  approximation  to  the  optimal  Bayes  procedure.  We  do  not  know 
if  this  procedure,  although  very  intuitive,  is  the  one  which  minimizes  the 
risk  for  the  squared  loss  function. 

Pinal  Remark: 


Note  that  the  criteria  used  depend  only  on  the  values  of  p^  and  p2, 
not  on  the  fact  that  we  have  an  exponential  distribution.  This  shows  the 
generality  of  the  methods  discussed.  On  the  other  hand,  suppose  that 
q2.  q^  »  1  -  (q^  ♦  q2)  are  the  probabilities  of  a  battery  failing 
respectively,  before  tj,  in  the  interval  (tj,  t2) ,  and  after  If  we 
represent  our  prior  opinion  about  (q^  q2<  q^)  by  a  Dirichlet  distribution, 
then  the  Bayes  solution  in  this  case  is  very  similar  to  the  one  presented 


